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ABSTRACT 

We construct a physically motivated model for predicting the properties of the rem- 
nants of gaseous galaxy mergers, given the properties of the progenitors and the orbit. 
The model is calibrated using a large suite of SPH merger simulations. It implements 
generalized energy conservation while accounting for dissipative energy losses and star 
formation. The dissipative effects are evaluated from the initial gas fractions and from 
the orbital parameters via an "impulse" parameter, which characterizes the strength 
of the encounter. Given the progenitor properties, the model predicts the remnant 
stellar mass, half-mass radius, and velocity dispersion to an accuracy of 25%. The 
model is valid for both major and minor mergers. We provide an explicit recipe for 
semi-analytic models of galaxy formation. 

Key words: galaxies: interactions - galaxies: evolution - galaxies: elliptical and 
lenticular, cD - galaxies: formation - methods: . 



1 INTRODUCTION 

Major mergers between galaxies ar e central to the formation 
, and evolution of ell iptical galaxies ("Toomr e fc Toomrd[l973 : 
' iToomrd [iOTtI : iMThos & Hornquist 1994). The hierarchical 
buildup of galaxies in the ACDM cosmology consists of a 
sequence of mergers, of which a significant fraction are "ma- 
jor," involving progenitors with a mass ratio larger than 1:3. 
The gravitational interactions in such mergers have a dra- 
matic effect on the dynamics and morphology of the galaxies, 
in particular turning rotating disks into pressure-supported 
spheroids. If the progenitors also contain gas, the mergers 
induce starbursts followed by gas consumption, which leads 
to aging stellar populations. The modeling of major merg- 
ers is therefore a key element in the attempts to confront 
the broad picture of galaxy formation with detailed obser- 
vations. This is commonly performed via simulations incor- 
porating Semi- Analytic Models (SAMs), where the complex 
physical processes are modeled using simplified parametric 
recipes. 

Advanced SAMs are currently attempting the non- 
trivial task of following the sizes and internal velocities of 
galaxies. For disk galaxies, sizes are evaluated using the halo 
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virial radii _Rvir and spin parameters A via -Rdisk — Ai?vir, 
with some modifications due to the halo density profile (Fall 
& Efstathiou 1980; Mo, Mao & White; Bullock, Dekel et 
al. 2001; Button et al. 2006). The sizes of the remnants of 
gas-poor ("dry") mergers, where the dominant interaction 
is gravitational, can be extracted from the properties of the 
progenitors and the orbital energy by assuming conse rvation 
of en ergy and relaxation to virial equilibrium (Cole et al.l 
|2000D . These considerations work well in simulations of rel- 
atively dry mergers, quite independently of the details of the 
orbit. 

While dry mergers m ay dominate in t he formation of 
the most massive galaxies (|Naab et al. ||200d), the most com- 
mon mergers are "wet" mergers of gaseous galaxies. It is 
thought that gas proc esses play an important role in the 
formation of ellipticals (^ Robertson et ah l l2006l : lDekel fc Coxl 
2006; Ciotti et al. 2007). As demonstrated below, the sizes 
predicted by dissipationless energy conservation can be off 
by a factor of a few for gas-rich mergers. Our goal is to con- 
struct a more accurate recipe to predict the size and velocity 
dispersion of the remnant of a wet merger given the prop- 
erties of the progenitors and the orbital parameters. Cos- 
mological mergers involve a complex mixture of variables 
(such as orbital parameters, gas fractions, mass ratio and 
bulge fraction) and physical processes (such as star forma- 
tion and feedback) , all of which can influence the properties 
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of the merger remnants. G iven a rich suite of high-reso lution 
SPH merger simulations l|Coxll2004l : ICox et al.lr2006h . that 
span the available parameter space and physical processes, 
albeit in a rather sparse and nonuniform manner, we seek 
a model that will properly represent the simulation results 
and enable an interpolation between them as well as an ex- 
trapolation to outside the simulated regime. 

For such a recipe to be successful, it should be based 
on a toy model that grasps the essence of the main physical 
processes involved in wet mergers. Our intuition is guided by 
the finding from the simulations that the remnants are more 
compact when the initial gas fraction is higher and when the 
first passage involves a stronger tidal impulse, namely when 
a larger fraction of the orbital energy turns into internal ki- 
netic energy. At a first glance, this may seem surprising, as 
a system that gains energy is not expected to become more 
tightly bound, and indeed, the remnants of dry mergers are 
not very sensitive to the strength of the impulse. This de- 
pendence on gas fraction and on the impulse implies that 
a higher gas fraction and a stronger impulse are associated 
with a higher degree of dissipation, via shocks, collisions of 
gas clouds, and induced gas flows toward the centres of the 
merging systems. The resultant higher gas densities enhance 
the energy losses to radiation, leaving behind a more tightly 
bound remnant. In parallel, the higher degree of dissipation 
yields a stronger burst of star formation, which tends to be 
focused in the central region of the remnant. This under- 
standing is the basis for our proposed recipe, which charac- 
terizes the merger by the impulse at first passage, evaluates 
the associated degree of dissipation and the resultant ra- 
diative energy losses and star formation, and accounts for 
these energy losses in the energy balance. A few free param- 
eters with values of order unity can hopefully compensate 
for the crude approximations made. These approximations 
include, for example, an assumption of structural homology 
between the progenitors and remnant. The physically moti- 
vated recipe is then calibrated using the merger simulations, 
and its success is to be judged by its accuracy in matching 
the simulated remnant properties. 

In 321we describe the simulations used for this study. In 
32]we present the details of our model for predicting remnant 
properties. In ^ we generalize the model to unequal mass 
mergers. In ^we discuss certain limitations of our model, 
and in !j6] we summarize our conclusions. Appendix A dis- 
cusses the details of our impulse approximation. Appendix 
B presents an explicit recipe for SAMs. 



2 MERGER SIMULATIONS 
2.1 Numerical Code 

The numerical simulations analyzed in this work are part of 
a large suite of galaxy merg er simulations designed to s tudy 
the induced star formati on (ICox 20041: ICox et al. 2006h and 
observable counterparts (| Jonsson l2006l ; 1jonsson et al.ll2006l ) 
of such events. Details of these simulations can be found in 
the above references, but we include here a brief description 
for completeness. 

All numerical simulations performed in this 
work use the N-Bo d y/SP H code GADGET 
l|Springel. Yoshida fc Whitel (iool). Hydrodynamics are 



included via the Lagrangian technique of smoothed particle 
hydrodynamics ( SPH). We use the "conserva tive entropy" 
version of SPH (|Springel fc Hernguistj |2003 ). Gas is as- 
sumed, for simplicity, to be a primordial plasma that can 
radiatively cool via atomic and free-free emission. 

All of the numerical simulations presented here include 
star formation. Stars are formed in regions of gas which are 
above a critical density for star formation at a rate propor- 
tional to the local gas density and inversely proportional to 
the local dynamical time-scale. The efficiency of star for- 
mation is fixed by requiring star formation to follow the 
observed correlat ion between gas and star-formation rate 
surface densities (|Kennicuttl ll998'). 

We also include a simple prescription to simulate the 
effects of feedback from massive stars. This feedback acts to 
pressurize the interstellar medium and regulates the conver- 
sion of gas to stars. Det ails of this model and the parameter 
choices can be found in lCox et al. I l|2006l ). Specifically, most 
simulations studied in this paper used the n2med parame- 
ter set. Under these assumptions the gas pressure increases 
as the density squared; i.e, star-forming gas has a "stiff" 
equation of state. Other cases are discussed in ^ 

The simulations presented here adopt a gravitational 
softening length h = 400 pc for the dark matter particles 
and 100 pc for the stellar and gas particles. We remind the 
reader that, in GADGET, forces between neighboring par- 
ticles become non-Newtonian for separations < 2.3 h. 

2.2 Initial Galaxies 

All of the simulations presented here are mergers between 
two identical disk galaxies, except for cases discussed in 
The disk galaxy models are motivated by observations of 
low-redshift galaxies. In some cases we made systematic 
studies of varying progenitor galaxy properties (e.g. vary- 
ing gas fraction in the G3 gas fraction series). While many 
of these varied cases would not look like typical low-redshift 
galaxies, we made no attempt to vary properties in such 
a way as to capture variation with redshift. Furthermore, 
the simulations are not cosmological since the two galax- 
ies are isolated. Disk galaxies are constructed in equilibrium 
and contain dark matter, an exponential stellar disk, an ex- 
tended exponential gas disk, and some contain a dense cen- 
tral bulge. Our suite consists of five main types of models, 
detailed in Table [1] 

(i) Z galaxies are gas-poor bulgeless disks and are roughly 
modeled after the Milky Way. 

(ii) D galaxies are of the mass of the Z's, are bulgeless 
disks, and have a high gas fraction. 

(iii) Y galaxies are of the mass of the Z's, are bulgeless 
disks, and have a high gas fraction. 

(iv) Sbc galaxies are modeled after local Sbc-type spirals, 
with a small bulge and high gas fraction. 

(v) G galaxies span a range of mass, bulge fraction, 
and gas fraction. Their properties are taken from statis- 
tical samples of local galaxies, including the SDSS. Their 
dark matter halos have not been adia batically contracted 
iBlumenthal et al.lll984lMo et all 19981 ). unlike all the other 
models. 

The ratio of the gas to stellar exponential radii varies 
with model type. For the Z, D, and Y models, the gas and 
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Table 1. Properties of progenitor galaxy models. Mtot is total 
mass, baryons plus dark matter; c is concentration (-Rvir/f's); 
Mstars is the initial stellar mass; B/D is the bulge-to-disk ra- 
tio; /g is the initial gas mass divided by Mtot; R1/2 initial 
three-dimensional stellar half mass radius. 



Type Mtot 


c 




B/D 


/g 


Rl/2 


(IQiOMo) 




(1O1OM0) 






(kpc) 


Milky Way Series 


D 1.4 


20 


0.036 





0.025 


1.16 


Y 14.0 


15 


0.3 





0.02 


2.85 


Z 143.0 


12 


5.1 





0.004 


4.04 


Sbc Series 


Sbc 81.4 


11 


4.92 


0.26 


0.066 


7.15 


G Series 


GO 5.0 


14 


0.1 


0.02 


0.012 


1.84 


Gl 20.0 


12 


0.5 


0.06 


0.010 


2.33 


G2 51.0 


9 


1.5 


0.11 


0.009 


2.90 


G3 116.0 


6 


5.0 


0.22 


0.011 


3.90 


G3 Gas Fraction Series 


G3gfl 116.0 


6 


3.6 


0.32 


0.023 


3.49 


G3gf2 116.0 


6 


2.6 


0.52 


0.031 


2.89 


G3gf3 116.0 


6 


1.5 


1.34 


0.040 


1.77 


G3gf4 116.0 


6 


5.3 


0.20 


0.005 


3.96 



stellar radii are equal. For the Sbc and G models, the gas 
radii are three t imes the st ellar radii. For more detail on 
these models see ICoxl (j2004l ). 

2.3 Merger Orbits 

The sizes of merger remnants are affected by the initial or- 
bits and orientations. To understand this relationship a suf- 
ficient exploration of the merger orbit and orientation pa- 
rameter space is required. To this end, we perform mergers 
on an identical orbit with various orientations of the merg- 
ing galaxies. We also perform mergers with many different 
orbits. Of all of the galaxy models, the Sbc models have the 
largest variety of orbits and orientations. The majority of 
the orbits in the suite are parabolic or near parabolic with 
eccentricities of 0.9 to 1.0. While these orbits are g enerally 
motivat ed by statistics from N -body simulations (|Bensonl 
[2005. : .Khochfar fc Burkertir2006i '). the distribution of orbits 
in the merger suite was not designed to exactly duplicate 
these statistics. 

3 MODELING PROPERTIES OF THE 
REMNANTS 

3.1 Non-dissipative Energy Conservation Model 

Energy conservation during the merger may be a useful con- 
straint to impose. If the progenitors and the remnant are 
homologous, then energy conservation and the virial theo- 
rem may be applied to the stellar systems as if they are 
self-gravitating without introducing a large error. Current 



SAMs emp loy such considerations to predict m erger rem- 
nant radii l|Cole et al.l I2OO0I : iHatton et all 120031 ). We start 
by summarizing the model of Cole et al. 

In order to take orbital energy into account, this model 
assumes that the baryonic components of the progenitors 
spiral in under dynamical friction, losing energy to the outer 
dark matter halo, until reaching a distance that equals the 
sum of their three-dimensional stellar half-mass radii. En- 
ergy conservation is assumed from this point on. Thus the 
orbital energy term in the conservation equation is equal to 
the energy of a zero eccentricity (circular) orbit of the two 
galaxies with a constant separation equal to the sum of their 
half-mass radii, 

{Ml + Ma)' _ M? Mj 1 M1M2 
Ri Ri R2 c Ri + R2' 

This equation relates the stellar three-dimensional half- 
mass radius of the remnant, Rf, to the masses and three- 
dimensional stellar half-mass radii of the progenitors. Mi, 
M2, Ri and R2, respectively. For major mergers, the masses 
are assumed to include the total stellar mass plus, as a rather 
arbitrary choice, twice the dark matter masses within Ri 
and R2 respectively. The constant c is a structural param- 
eter which relates GM /R to the actual internal binding 
energy. It is assumed to take the same value, c — 0.5, for 
both the progenitors and the remnant. 

This model, though crude, provides a general framework 
for estimating the outcome of a given merger. However, it 
has not been previously tested against realistic simulations 
that include gas dynamics and star formation. We begin by 
applying this recipe to the cases in our merger simulation 
suite. The predictions of the model are plotted against the 
actual half-mass radii of the simulated remnants in Figure [T] 
Throughout this paper we use the fractional rms scatter, 5", 
to assess goodness of fit. 

r, /I (fprcdictcd ^ ftruc)' 

^=v^^ — p^'^ — ' 

where P is some property that we are trying to predict from 
the initial conditions and are measuring from the simulation 
for comparison. 

The predictions of this simple model deviate from the 
true radii by 5* = 0.50. Figure [1] shows that the model sys- 
tematically over-predicts the radii of the remnants. This is a 
straightforward result of ignoring the radiative energy losses. 
Note that the predictions are best for the cases where the 
progenitors have the lowest gas fraction, cases Z and G3. 
For each type of progenitor, we see a wide spread in actual 
sizes for a given predicted size. This results from the varia- 
tions in orbits and orientations, which are not addressed by 
the model. Despite the shortcomings of the dissipationless 
energy conservation model, it is a useful starting point. Our 
next goal is to correct for the dissipative effects. 

3.2 A Toy Model for Radiative Losses 

While the above dissipationless model works quite well in the 
case of low gas fraction, dissipative losses are likely to play an 
important role in the gas-rich mergers that were especially 
frequent in the early epochs of galaxy formation. In particu- 
lar, they seem to be a crucial element in the formation of the 



4 Covington et al. 



Table 2. The orbital parameters of each of the simulations used in this study. Values are calculated from the orbital initial conditions 
assuming a point mass orbit. Rperi denotes pericentric distance, e denotes eccentricity. For equal mass mergers 8i and 02 denote the 
orientation of the first and second galaxy with respect to the orbital plane, where for 9i = the first galaxy is aligned with the orbital 
plane. For unequal mass mergers 9 denotes the orientation of the orbit of the smaller progenitor with respect to the orientation of the 
larger progenitor. Again, 6 = represents a case where the merger orbit aligns with the larger progenitor. 



Major Mergers (55 runs) 



name 




e 


fi 


02 


name 




Wl 


0^ 
02 


Dlmf-ul 


1.1 




1.0 





30 


Sbc201-u4 


11.0 


1.0 





30 


D2mf-ul 


2.4 




1.0 





30 


Sbc202-u4 


11.0 


1.0 


180 


30 


D3mf-ul 


4.1 




1.0 





30 


Sbc203-u4 


11.0 


1.0 


180 


210 


D4mf-ul 


6.4 




1.0 





30 


Sbc204-u4 


5.5 


1.0 





30 


D5mf-ul 


8.9 




1.0 





30 


Sbc205-u4 


44.0 


1.0 





30 


D6mf-ul 


12.9 




1.0 





30 


Sbc206-u4 


11.0 


1.0 


90 


30 


D7mf-ul 


0.6 




1.0 





30 


Sbc207-u4 


11.0 


1.0 


270 


30 


GOGOa-ul 


2.24 




0.95 


-30 


30 


Sbc208-u4 


5.5 


1.0 


180 


30 


GOGOr-ul 


2.24 




0.95 


150 


30 


Sbc209-u4 


5.5 


1.0 


180 


210 


GOGO-ul 


2.24 




1.0 


-30 


30 


Sbc211-u4 


44.0 


1.0 


180 


210 


GlGla-ul 


2.96 




0.95 


-30 


30 


Sbc212-u4 


11.0 


0.9 





30 


GlGlr-ul 


2.96 




0.95 


150 


30 


Sbc213-u4 


25.0 


0.8 





30 


G2G2r-ul 


3.82 




0.95 


150 


30 


Sbc214-u4 


44.0 


0.8 





30 


G2G2-U1 


3.82 




0.95 


-30 


30 


Sbc215-u4 


100.0 


1.0 





30 


G3blv5G3blv5-ul 


13.6 




0.95 


-30 


30 


Sbc216-u4 


100.0 


0.8 





30 


G3G3a-ul 


13.6 




0.95 


-30 


30 


Sbc217-u4 


11.0 


1.0 


90 


90 


G3G3b-ul 


13.6 




0.95 


-30 


30 


Sbc218-u4 


11.0 


0.9 


180 


210 


G3G3r-ul 


13.6 




0.95 


150 


30 


Z2m-ul 


7.1 


1.0 





30 


G3gflG3gflb-ul 


13.6 




0.95 


-30 


30 


Z7m-ul 


21.4 


0.9 





30 


G3gf2G3gf2b-ul 


13.6 




0.95 


-30 


30 


Z8m-ul 


35.7 


0.8 





30 


G3gf3G3gf3b-ul 


13.6 




0.95 


-30 


30 


Z9m-ul 


1.7 


1.0 





30 


G3gf4G3gf4b-ul 


13.6 




0.95 


-30 


30 


ZlOm-ul 


3.9 


1.0 





30 


Ylmf-ul 


2.9 




1.0 





30 


Zllm-ul 


14.2 


1.0 





30 


Y2mf-ul 


5.7 




1.0 





30 


Z12m-ul 


22.2 


1.0 





30 


Y3mf-ul 


10.0 




1.0 





30 


Z13m-ul 


30.4 


1.0 





30 


Y4mf-ul 


15.7 




1.0 





30 


Z14m-ul 


44.7 


1.0 





30 


Y5mf-ul 


21.4 




1.0 





30 


Z15m-ul 


8.1 


1.0 





30 


Y6mf-ul 


31.4 




1.0 





30 












Minor Mergers (22 runs 


name Rperi (kpc) 


e 


e 


Mass Ratio 




name Rj 




e 





Mass Ratio 


G1G0-U3 2.96 


0.95 


-30 




1:3.9 




G3Glf-ul 


27.2 


0.95 


-30 


1:5.8 


GlG0r-u3 2.96 


0.95 


150 




1:3.9 




G3Glg-ul 


64.4 


0.95 


-30 


1:5.8 


G2G0-U3 3.82 


0.95 


-30 




1:10 




G3Glh-ul 


120 





-30 


1:5.8 


G2G0r-u3 3.82 


0.95 


150 




1:10 




G3G2a-ul 


13.6 


0.95 





1:2.3 


G2G1-U3 3.82 


0.95 


-30 




1:2.6 




G3G2b-ul 


13.6 


0.95 


-90 


1:2.3 


G2Glr-u3 3.82 


0.95 


150 




1:2.6 




G3G2c-ul 


13.6 


0.95 


-60 


1:2.3 


G3Gla-ul 13.6 


0.95 







1:5.8 




G3G2d-ul 


13.6 


0.95 


180 


1:2.3 


G3Glb-ul 13.6 


0.95 


-90 




1:5.8 




G3G2e-ul 


6.8 


0.95 


-30 


1:2.3 


G3Glc-ul 13.6 


0.95 


-60 




1:5.8 




G3G2f-ul 


27.2 


0.95 


-30 


1:2.3 


G3Gld-ul 13.6 


0.95 


180 




1:5.8 




G3G2g-ul 


64.4 


0.95 


-30 


1:2.3 


G3Gle-ul 6.8 


0.95 


-30 




1:5.8 




G3G2h-ul 


120 





-30 


1:2.3 



Fundamental Plane of elliptical galaxies (jRobertson et all 
[2OO6; Dokcl & Cox 200QV We therefore wish to incorporate 
the radiative losses in our model. 

When two gas-rich galaxies merge, a number of pro- 
cesses cause gas interactions and result in radiative energy 
losses. Tidal torques during a close pass can decrease the 
angular momentum in the gas disks and induce inflows into 
the galaxy centres. In a nearly radial encounter, the gas disk 
of one galaxy collides with that of the other creating shocks, 
which can result in loss of angular momentum. Furthermore, 



tidal forces during a merger introduce orbital crossings and 
density perturbations within the gas disks which ultimately 
lead to an increased gas collision rate. As the gas clouds 
collide and radiate away their kinetic energies, they fall to- 
ward the centre of the galactic potential well. This results 
in higher gas densities which lead to star formation and fur- 
ther radiative losses. Hence, the orbit, energy losses, and 
star formation are intimately linked. 

The first step toward predicting energy losses from a 
given initial configuration is to characterize the perturbative 
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Figure 1. Dissipationless energy conservation model. Tlie tliree- 
dimensional lialf-mass radii predicted by the model versus the 
actual radii in the simulated remnants. The symbols and colors 
denote the type of progenitors (Table 1): D (black crosses), G 
(purple asterisks), Y (blue diamonds), Z (green triangles) and Sbc 
(red squares). The symbols for G's have four different sizes which 
represent the different types of G progenitors. Larger symbols 
represent more massive galaxies. With a relatively large scatter 
of 5 = 0.5, this dissipationless model systematically over-predicts 
the sizes of the remnants, because it does not address the radiative 
energy losses. 



strength of a given encounter. More specifically, we need to 
assess the dependence of the gas collisions and dissipation 
on the properties of the orbit. We find that the dissipation 
resulting from a given orbit can be characterized using the 
orbital parameters at the first close pass. This may seem 
surprising at first, since often most of the stars are formed 
in the final coalescence rather than at the initial encounter, 
but it can be explained by the sequence of events following 
the first pass. The gas disks are perturbed during the first 
pass, inducing a continuous gas in-fall toward the galactic 
centre, lasting ~ 1 Gyr. A stronger disturbance in the first 
pass leads to a larger buildup of gas in the progenitor cen- 
tres. Some of this gas is involved in a first starburst imme- 
diately following the first pass, but a large fraction of this 
gas serves as a reservoir for star formation during the later 
stages of the merger, especially the violent final coalescence. 
Consequently, an orbit that is more disruptive on the first 
pass also suffers more energy losses, and forms more stars, 
during the final coalescence. This is illustrated in Figure [2] 
During a close pass of two galaxies, orbital energy is in- 
jected into the internal kinetic energies of both galaxies. We 
define this "impulse" as the difference between the peak in 
the total internal kinetic energies during the encounter and 
the total initial internal kinetic energies of the two progeni- 
tors. We approximate the impulse on galaxy 1 by 



AE : 



AGHdl,,,M2. 



1/2 

pen 



i(^pcri ~^ ^ -Rl,tot -Rpcri + C R-^^^ 



(3) 



where Mi^tot and M2,tot are the total mass of the perturbed 
and perturbing galaxies respectively, baryons plus dark mat- 
ter, and i?i,tot is the total half mass radius of the perturbed 
galaxy. 7?peri is the "pericentric distance" of the first passage, 
as calculated from the initial orbit assuming point masses. 



Rp„i = 5.5 kpc 
Rp„i = 1 1 kpc 
R„„i = 100 kpc 



o 0.00 [ 
0.0 



.0 1.5 

Time (Gyr) 



Figure 2. The central gas density in the inner 0.5 kpc during 
three separate Sbc mergers with initial orbital pericentric dis- 
tances of 5.5, f 1 and 100 kpc (solid black, dashed blue, and dash- 
dot red respectively). In all cases, the initial pass occurs at 0.6 
Gyr. Final coalescences occur at roughly 1.5, 1.8 and 2.4 Gyr re- 
spectively. The plot demonstrates that the central gas density in 
the period following the first pass is increasing with decreasing 
pcricentre distance, and that the density at the final coalescence 
is following the same trend. 




10" 1 0^ 

True Impulse 

Figure 3. The approximation for the impulse against the mea- 
sured impulse from the simulations. S = 0.40. 



The best-fitting values of the parameters are found to be 
A = 1.6, B = 1.0 and C = 0.006. The fit of the equation 
to the impulse measured from the simulations in shown in 
Figure [3] A more detailed discussion of this impulse approx- 
imation is given in Appendix A. 

We characterize the dissipative strength of a merger us- 
ing the fractional impulse at the first pass, fk = AE/Ktot, 
where AE is the impulse and Ktot is the total initial inter- 
nal kinetic energy of the galaxy, baryons and dark matter. 
We find that /new, the fraction of new stars formed in the 
merger relative to A/tot, is indeed proportional to the frac- 
tional impulse and to the gas fraction /g. 



/n 



C'ncw /g /k. 



The best fit to our simulations is obtained for Cnc 



(4) 
0.3. 
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Figure 4. Fraction of new stars formed during the merger: the ap- 
proximation based on gas fraction and impulse versus the actual 
value in the simulations. A stronger perturbation from a larger 
impulse enhances the collision rate between gas clouds, and thus 
induces more star formation. 



1.00 




Figure 5. Fractional energy loss. The approximation of radiative 
losses in equation[6l which is based on gas fraction and impulse, is 
shown against the actual radiative energy losses. It demonstrates 
a crude proportionality, with the line plotted y = OAx. 



The gas fraction determines the overaU norrnahzation of star 
formation, whereas the impulse factor tracks the variations 
with orbit. This prediction for the fraction of new stars is 
plotted in Figure |4] against the actual fraction of new stars 
in the simulated remnants. In order to consistently treat 
the star formation in simulations that were run for dif- 
ferent lengths of time, and in keeping with previous work 
l|Somerville et al1l200ll ). we separate star formation into a 
merger-induced burst component and a quiescent compo- 
nent. We accomplish this by subtracting out star formation 
measured from simulations of each progenitor in isolation. 
The quantity /now in equation |3] refers to the burst compo- 
nent only. 

Given the above correlation between impulse and star 
formation, we proceed to predict the energy loss by assuming 
that the colliding gas that forms new stars during the merger 
loses a constant fraction of its kinetic energy in the collision. 
In more detail, we make the following assumptions: 

(i) Gas clouds have initial velocities approximately equal 
to the average initial velocities of the dark matter/baryon 
system. 

(ii) The average impulse per mass imparted to gas that 
will form new stars is approximately equal to the average 
impulse per mass for the entire dark matter/baryon system. 

(iii) The energy lost during the merger is proportional to 
the kinetic energy of the gas that will form stars. This in- 
cludes: 1) the initial kinetic energy, and 2) the kinetic energy 
gained from the impulse. 

While simplistic, these assumptions allow us to make 
a connection between radiative losses and star formation. 
Assumption (i) would be valid if the system resembled an 
isothermal sphere. Assumption (iii) provides a sensible way 
for parameterizing the energy loss, given that it is associated 
with collisions of gas clouds. Under this assumption, the 
radiative energy loss can be written as 



The first term is the total initial internal kinetic energy of 
the gas that forms stars, and the second term is the energy 
imparted to that gas by the impulse. Dividing by Ktot, and 
using equation |4] gives 



/rad OC /ncw(l + /k) OC /g /k (1 + /k). 



(6) 



This expression is plotted in Figure [S] against the actual 
radiative energy losses in our simulations, showing a crude 
proportionality, with the proportionality constant ~ 0.4. 



3.3 A Dissipative Model for Remnant Radii 

Guided by this toy model for dissipative losses during a 
galaxy merger, we construct a modified energy conservation 
equation for predicting remnant radii. Note that the final 
stellar mass can be computed using the initial mass and the 
equation for star formation @. We choose to include in our 
energy formula the initial stellar mass plus the mass of the 
stars formed during the merger, such that the dissipative 
term tracks the loss of energy of the progenitor's gas that 
becomes stars in the remnant. The energy equation is 

-Eint.f = -Eiiit.i + -Brad + -Eorb , (7) 

where the final and initial internal energies are 

CintG(A/l + Ma + Mnow.l + A/new,2)^ 

-t/int.f ~ t: (,oj 



EiYit,! — — CintG 



(Ml + Mnc^.i)'-^ , (Ma + Mnow.a) 



+ 



R2 



(9) 



-Erad OC (/new A'tot + /new AE) . 



(5) 



and Ml, Ma, Ri, and R2 are the initial stellar masses and 
the corresponding three-dimensional half-mass radii of the 
progenitors, Rf is the final three-dimensional half-mass ra- 
dius of the remnant, and Mnow.i and Mnow.a are the mass of 
new stars formed during the merger in galaxies one and two 
respectively, as predicted by our model. The constant struc- 
tural parameter, which relates the internal energy of the 
systems to GA4^ /R, is determined by best fit to the radii of 
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the simulated remnants to be Cint — 0.5. In assuming ho- 
mology, we set the structural constants for the progenitors 
and remnants to the same value. 

The radiated energy term is the sum of the losses from 
the two progenitors, 



-CradSi^lA'i/g,i/k,i(l + /k 



(10) 



where Ki, fg^i and /k,i are the initial internal kinetic energy, 
gas fraction and impulse corresponding to progenitor i. The 
constant of proportionality relating our expression for en- 
ergy loss to the actual energy lost is determined by best fit 
to the radii of the simulated remnants to be Crad — 1-0. 
The orbital energy is: 

G(Ml + Mnow,l)(M2 + Mncw,2) 



Eorb — 



+ -(Ml + Mncw,l)Vi^ + -(M2 + Mnc«,2)Vi 

where i?scp, the distance between the progenitors' centres of 
mass, and Vi and V2, the centre of mass velocites of each 
progenitor, are defined at the beginning of the encounter. 
Most of our simulations are nearly parabolic, such that this 
term is close to zero. 

Our model predictions for the remnant radii are plot- 
ted against the actual stellar half-mass radii of the simulated 
merger remnants in Figure El The overall scatter is 5 = 0.21, 
which is a significant improvement over the dissipationless 
model. There is no obvious systematic error, and the fit is 
good both for the wide range of progenitor properties and 
for the different choices of orbital parameters. The spread 
for a particular progenitor type results from differences in 
orbit and orientation. For most cases this orbital spread is 
well fit by the model. This is particularly evident with the 
D and Y series, for which only orbit is varied and not ori- 
entation. The spread in Sbc series is not quite as closely 
tracked by the model. This additional spread is the result of 
differing orientations. Orientation is not taken into account 
in the model for reasons described in i]5.2l We note that the 
performance of the final model for radii is somewhat better 
than what might have been expected from the quality of the 
prediction of radiative losses, shown in (5] However, this is 
not as surprising if one remembers that radiative losses are 
only a correction on the non-dissipative model. Even if this 
correction is not exact, it always acts in the proper direction 
and makes the remnants more compact. 

3.4 Velocity Dispersion 

iNaab et al.l l|2006t ) have shown that the kinematics of merger 
remnants change as a function of initial gas fraction. When 
predicting the central velocity dispersion of merger remnants 
we take gas fraction into account in two ways: first by using 
the gas fraction dependent predicted radius, and second by 
adjusting the central dark matter fraction to account for 
the rearrangement of gas and subsequent formation of new 
stars. In order to compute the stellar velocity dispersion of 
the merger remnants, we implement a virial relation of the 
type 

GMdyn 



Cvir 



R 



(11) 



where Mdyn is the dynamical mass of the system. This mass 
includes all of the stellar mass of the system and also a 



contribution from the dark matter near the centre of the 
galaxy. Cvir is a constant that varies slightly with galaxy 
structure and also accounts for the conversion between the 
three-dimensional radius and the projected, two-dimensional 
velocity dispersion. 71 is a characteristic radius (e.g., the stel- 
lar half-mass radius), cr is a line-of-sight velocity dispersion 
of stars inside the projected half-mass radius. We measure 
this from the simulations by averaging over 50 random pro- 
jections. 

Given the predicted values for R{ and the final stel- 
lar mass, one may attempt to estimate a^. However, varia- 
tions in central progenitor dark matter fractions and differ- 
ing amounts of star formation result in significantly different 
ratios of Mdyn to the total stellar mass in the remnants. Be- 
cause we are dealing with an approximate virial relation of a 
system that actually is not completely self-gravitating, it is 
not obvious from first principles what the exact contribution 
to Afdyn from dark matter should be. We track the dark mat- 
ter contribution by estimating the dark-matter fraction for 
each remnant. Therefore the uncertainty above, concerning 
which dark matter mass to include, translates into an uncer- 
tainty in which radius to choose for defining a dark-matter 
fraction. Because we use the three-dimensional stellar half- 
mass radius in our virial relation, it would be reasonable 
to choose a radius ^ Rf. However, one can pick a range of 
radii and still achieve sensible results by making slight ad- 
justments to Cvir. We find the best results when we focus on 
the dark matter fraction inside half of the three-dimensional 
stellar half-mass radius. We define the dark-matter fraction 
within a given radius by 



(12) 



(Mdm + Matars)' 

where Mdm and Mstars are the dark matter and stellar 
masses inside that radius, respectively. Much of the varia- 
tion in the remnant dark matter fraction, /dm,f , is due to the 
variation in initial dark-matter concentrations and baryon 
distributions. However, another important effect is the ten- 
dency of new stars to form near the galaxy centre, causing 
mergers that form more stars to end up with lower /dm,f 
values. We thus predict the final dark-matter fraction using 
the initial dark matter masses and our model prediction for 
the mass in new stars, 

m,2 



/dm,f = 



(13) 



A/dm,l + A/dm,2 + atars(Afl + Mi + A4"new ) ' 

Mdm.i and A/dm, 2 are the dark matter masses inside half 
of the three-dimensional stellar half-mass radii of the pro- 
genitors. All and A/2 are the stellar masses of the progeni- 
tors, and A/ncw is the total mass of stars formed during the 
merger as predicted by equation This expression simply 
assumes that the inner region of the remnant contains the 
same amount of dark matter as the sum of the inner regions 
of the progenitors, and that a fixed fraction, Cstars, of the 
final stellar mass is inside one-half of the three-dimensional 
stellar half-mass radius. We find that the best fit to the sim- 
ulated remnants is Cstars — 0.35. 

Our modified virial relation becomes 

G(Afi + Ah + A/„ow) 



Cvir 



(14) 



-Rf(l - /dm.f) 

The best-fitting value from our simulations is Cvir — 0.30. 
The model predictions for the stellar line-of-sight velocity 
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dispersions are shown in Figure [6l in comparison with the 
simulated values. The scatter in cr^ is 5 = 0.24, and the 
overall performance is similar to that of the model for the 
final radii. There are no obvious systematic errors, and the 
predictions properly capture the variations due to either pro- 
genitor properties or orbital parameters. 



4 UNEQUAL MASS MERGERS 

The cases studied so far were mergers of identical galaxies, 
and the model was built using the simulations of these cases. 
However, most mergers in the universe are of unequal mass 
galaxies. More recent sim ulations were run of unequal mass 
mergers l|Cox et al.ll2007f ). and given the importance of un- 
equal mass mergers in the real universe we here examine the 
accuracy of our model for different progenitor mass ratios. 
For this study we used simulation s of unequal mas s merg- 
ers from the G series of galaxies l|Cox et al.ll2007^ . These 
mergers cover a range of mass ratios and a variety of orbital 
parameters. For unequal mass mergers, the predictions of 
remnant properties remain relatively accurate with 5* = 0.13 
about the previously fit relation for radii, and S = 0.20 for 
(see Figure (Tjl. Note that in some cases the spread does 
not fall directly along the line. This is the result of a number 
of simulations with the same progenitors and orbital param- 
eters, but with different orientations of the progenitors with 
respect to the orbital plane. Variations due to orientation 
do not have clear enough systematics to be accounted for in 
the model. We discuss this limitation in ii5.2l 

In major mergers, which have mass ratios greater than 
1:3, the gas disks of the galaxies are severely disrupted, and 
we see large flows of gas toward the galactic centres. How- 
ever, in minor mergers, which have mass ratios less than 1:3, 
these flows are much less pronounced because the gas disk of 
the larger progenitor is only modestly disrupted. For exam- 
ple, central gas densities in a major merger of two G3 galax- 
ies reach values roughly four times that of a minor merger 
(1:6) of G3 and Gl galaxies on the same orbit. Hence, dis- 
sipation is expected to be much less important in the cen- 



tres of the big progenitors of minor mergers. Our model 
captures this effect through the impulse dependence of the 
radiation term. For minor mergers where Mbig 3> A/smaii, 
the radiation term becomes insigniflcant in comparison to 
the internal energy term. It is also worth noting that the 
the star formation equation @ under predicts star forma- 
tion in the smaller progenitor in this extreme minor merger 
regime. However, this has little effect on the predictions of 
remnant properties since the mass of new stars formed in a 
much smaller progenitor is much smaller than the mass of 
the larger progenitor. 



5 CAVEATS 
5.1 Feedback 

Arguably the largest uncertainty in the physics of our simu- 
lations is in the prescription for feedback. The merger sim- 
ulations fit by our model all use the same prescription for 
feedback with the same parameters, hence, it would be use- 
ful to vary these feedback parameters and examine the effect 
on our model. The feedback recipe is characterized by two 
parameters: feedback efficiency and an equation of state pa- 
rameter, n, which sets the polytropic index. The effective 
pressure is 

P.^-p^+f"/^). (15) 

We explore three equations of state, n = 0, n = 1, and 
n = 2, where n = corresponds to an isothermal equa- 
tion of state and n = 2 results in a stiff equation of state. 
The feedback efficiency determines how quickly the feed- 
back energy is allowed to thermalize. Higher efficiencies re- 
sult in a quicker dissipation of the feedback energy. For each 
value of n, we examine two different values of efficiency. 
As a lower limit, we use a low efficiency value that gives 
just enough pressure to stabilize the disk. We also simulate 
cases of super-stable disks where the efficiency is set to ten 
times that needed for disk stabilization. In table 2, labels 
of the feedback parameter sets denote the values of n and 
efficiency chosen. Unstable cases with even lower feedback 



Predicting the Properties of Merger Remnants 9 




Table 3. Effect of differing feedback parameters on remnant 
properties. All simulations in the table are variations of a fidu- 
cial Sbc merger. Afnowstars is the total mass in new stars formed 
during the 3 Gyr simulation, burst and quiescent. R{^i/2 is the 
stellar three-dimensional half-mass radius of the remnant. Cf is 
the velocity dispersion of the remnants measured inside of the 
projected half-mass radius averaged over 50 random projections. 



Model 


-^^ncwstars 








(IQlOM©) 


(kpc) 


(km s~-^) 


nOmed 


5.75 


5.37 


253 


nOhigh 


2.68 


5.72 


218 


nlmed 


5.30 


5.34 


252 


nlhigh 


2.57 


5.75 


190 


n2med(fiducial) 


5.03 


5.62 


225 


n2high 


1.61 


6.95 


178 



efficiency were run but are not discussed here. Thus, for con- 
sistency our "low" efficiency cases are labeled "med"(e.g. 
n2med means n = 2 and the feedback efficiency is low). 
In general, with lower feedback, the maximum star forma- 
tion rate is higher, but the starburst duration is shorter. 
For a more detailed d escription of our feedback model see 
ICox et al.1 (|200i,|2003). 

To examine the effects of feedback on our remnant 
model, we take the fiducial Sbc merger and simulate identi- 
cal initial conditions with varied values of the feedback pa- 
rameters. As expected, feedback has a significant effect on 
star formation. The effect is twofold. First, both increased 
feedback efficiency and, to a smaller extent, stiffer equations 
of state (higher n) decrease the total number of stars formed. 
This can be a dramatic effect. As shown in Table |3] the to- 
tal number of stars formed varies by more than a factor of 
three. Secondly, feedback can change the radial distribution 
of stars formed. A stiffer equation of state will lea d to more 
stars being formed at large radii (|Cox et al.ll2006l ). 

Both of these effects can act in concert to increase the 
size of the merger remnant. However, the effect is weak 



enough that it takes both high feedback efficiency and a 
stiff equation of state to significantly alter the final half- 
mass radius (see Figure |8]). In our set of feedback models, 
the only model for which the radius is significantly different 
than the others is n2high. Even in this model the remnant 
radius is only a factor of ~ 1.3 larger than the most com- 
pact remnant. Differences in radius and compactness will 
also affect the central velocity dispersion. Higher efficiency 
and stiffer feedback models produce somewhat lower a. This 
effect is similar in magnitude to the effect on radius, with 
a factor of ~ 1.4 between the largest and smallest a. For 
velocity dispersion, our fiducial model lies near the centre of 
the distribution (see Figure|8]). These differences are roughly 
comparable to the scatter in our remnant model. 

While we have chosen a specific feedback model as a 
fiducial, and used it to run most of our cases, we note that 
there are currently no theoretical or observational motiva- 
tions for this choice. Feedback in the real universe could 
resemble any of these models. If feedback in the real uni- 
verse is significantly different from our fiducial case, then it 
is possible that the parameters of our model would require 
some tuning. Cnew is the parameter most affected by differ- 
ing feedback. The variations in star formation suggest that 
Cnow could vary roughly from 0.1 to 0.35. It is also possible 
that the structural parameters Cint, Cvir, and Cstars would 
need to be adjusted. However, aside from the mass of new 
stars, the relatively drastic changes in feedback produce rel- 
atively minor changes in remnant properties. This suggests 
that any tuning needed to match the real universe would be 
small. Further verification of the model will require testing 
with independent feedback recipes. 

5.2 Orbits 

There are several potential sources of uncertainty concern- 
ing our treatment of the initial orbits of the progenitors. The 
first regards the orbital ene r gy te rm in our energy conser- 
vation equation. ICole et al] (|200G| ) assume in their dissipa- 
tionless model that some of the orbital energy is transferred 
to the outer dark matter halo through dynamical friction. 
In order to estimate the energy transferred from the orbit 
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Figure 8. Various feedback models. The predicted values of radius (a) and velocity dispersion (b) for all of the variations on the 
feedback model. Large symbols represent high feedback efficiency, small symbols represent low efficiency, and squares, crosses, and 
triangles represent n=2, n=l, and n=0 respectively. 



to the baryonic remnant, they assume that the two galaxies 
come into a circular orbit at a radius that equals the sum 
of their half-mass radii. However, the orbits in our simula- 
tions actually tend to quickly radialize after the first close 
passage. We therefore adopt the better approximation that 
all the orbital energy ends up in the baryonic remnant. This 
too may be an over-simplification, as the energy deposited 
in the dark halo is likely to depend on the initial eccentricity 
of the orbit, and may get larger for increasingly tangential 
orbits. Unfortunately, our suite of simulations is not partic- 
ularly well-suited to exploring this dependence, as most of 
them start with initial orbital eccentricities close to unity. 
Our suite does contain three Sbc simulations and one Z sim- 
ulation with e = 0.8, and also one unequal mass merger, a 
G3-G2, with an initially circular orbit, e = 0. We find that 
these few cases are well fit by our model as is. However, given 
our limited sampling of non-radial orbits, we keep in mind 
that our orbital-energy term may require some modification. 
On the other hand, we note that statistical studies of cos- 
mologi cal merger orb i ts suggest that they tend t o be rather 
radial l|Bensoni [2OO5I : iKhochfar fc Burker^ 120061 ). and that 
all attempts we made to include a term for energy transfer 
to the dark matter did not significantly improve the model. 

A second potential concern is the initial orbital parame- 
ters we use. Our simulations typically start when the centre 
of one galaxy is near half of the virial radius of the other 
galaxy, whereas the typical analysis of halo or bits in N-body 
simulations refer to orbits at the virial radius (|Bensonll2005l : 
iKhochfar fc Burkertl200^ '). Furthermore, our simulations in- 
volve the final entry into this radius, whereas in cosmological 
situations many of the halos enter and leave the virial ra- 
dius more than once before merging. One practical concern 
is how to convert between orbits at 0.5-Rvir and at i?vir. 
Dynamical friction undoubtedly induces some evolution in 
the orbital parameters as the galaxies fall between 0.5 -Rvir 
and -Rvir- This effect should be evaluated using cosmological 
N-body simulations. 

Finally, our model does not take into account the effect 
of the initial orientation of the progenitors with respect to 
the orbital plane. Our suite of simulations does include a 



number of different initial orientations, including both pro- 
grade and retrograde encounters. Previous authors have ob- 
served systematic differences between encounters with dif- 
ferent orientations (|Naab et al.ll2006l : lBarnesll2002l '). but the 
effects on the remnants in our simulations do not show clear 
enough systematics for us to effectively characterize them. 
The performance of the model could be improved by in- 
cluding a term that depends on orientation, but the modest 
improvement did not seem to justify the inclusion of yet an- 
other parameter that would increase both the complexity of 
the model and the risk that the model is over- fit. 

5.3 Gas Fraction 

Our merger suite primarily consists of disky progenitors with 
high gas fractions. It is possible that we are missing some 
of the effects that domi nate the results of dry mergers of 
spheroids. For example, iBovlan-Kolchin et al.l (i2006) find 
that more radial encounters of dry ellipticals result in larger 
remnants. However, we see no evidence for such a trend, 
not even in our lowest gas fraction simulations, the Z series, 
where Mgas/ Mtaryon ~ 0.1. In fact, for these low-gas cases, 
we see essentially no systematic trend of remnant size with 
orbit, and a generally smaller scatter in remnant size. We 
conclude that if we miss any effect of this sort, it is likely 
to be small. Apparently, the dissipationless energy conserva- 
tion equation becomes a relatively accurate approximation 
in the limit of low gas fraction, and our model converges to 
the dissipationless model when the progenitors are gas-free. 

Another limitation of our model is that it is unlikely to 
apply to extremely gas-rich mergers, where Mgas ^ M stars ■ 
For such cases the initial energy term should account for 
the size of the gas disk, whereas our model only takes into 
account the initial size of the stellar disks. The highest gas 
fraction case that we include has a ratio of Mgas/Matars = 
3.1. For the cases that we simulate, where Mgas ^ Afatars 
we find that the final sizes are relatively insensitive to the 
initial sizes of the gas disks. Specifically, our suite includes 
cases where the gas disk scale length is equal to the stellar 
disk scale length and cases where the gas disk scale length is 
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three times the steUar disk scale length. Both of these cases 
are well fit by our model. 

One final note is that our merger simulation suite does 
not attempt to model galaxies at higher redshift by system- 
atically changing gas fraction and concentrat ion as has been 
done in other studies l|Robertson et al.ll2006f ). Since progen- 
itors are primarily modeled after low-redshift galaxies, this 
is the regime where the model is most well-calibrated. How- 
ever, the model is based on simple physical principles and 
is robust to systematic changes in gas and bulge fraction. 
We see no obvious reason why it should not apply to binary 
mergers at higher redshifts. 



6 CONCLUSIONS 

We have developed a simple toy model for the physical pro- 
cesses involved in wet mergers of galaxies, and have cali- 
brated it using a suite of hydrodynamical merger simula- 
tions. This modeling helped us to gain a better understand- 
ing of these processes, and provides a practical semi-analytic 
recipe for predicting post-merger galaxy properti es in SAMs. 

Cr ude models of t his sor t have been used b y Cole et al] 

l|2000l ). iHatton et al] l|2003f) . and IShen et al.l l|2003l ), but 
these models did not account for energy losses through dissi- 
pative processes, and they have not been calibrated against 
realistic merger simulations. Using a suite of merger sim- 
ulations, we have demonstrated the key role of dissipative 
energy losses in determining the final radii and velocity dis- 
persions of merger remnants. We found that the dissipative 
effects depend on the initial orbits of the progenitors. More 
violent, lower angular momentum orbits create greater dis- 
turbances in the gas disks, which in turn radiate more energy 
and produce more stars. This orbital "violence" can be pa- 
rameterized through an impulse approximation for energy 
exchange between the orbital and internal components dur- 
ing the first close pass of the encounter. 

We present a physically-motivated, simulation- 
calibrated model that is capable of predicting star 
formation, central dark matter fraction, remnant radius 
and remnant velocity dispersion, given the properties of the 
progenitors and the initial orbital parameters of a merger. 

The non-dissipative energy conservation model often 
predicts radii that are off by a factor of ~ 2 — 3, and it 
does not reproduce the spread due to orbital variations. Our 
model, which accounts for the dissipative energy losses, re- 
sults in only ~ 25% errors in the predicted radius and ve- 
locity dispersion when a wide variety of progenitor types is 
considered. For a given progenitor type, the error in remnant 
properties is reduced to ~ 10%, indicating that our model 
correctly captures the variation of remnant properties due 
to merger orbit. 

Since we used the whole available simulation suite to 
calibrate our model, via a few proportionality constants of 
order unity, a proper evaluation of the model performance is 
yet to be pursued using an independent suite of simulations. 
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7 APPENDIX A: IMPULSE APPROXIMATION 

A number of previous researchers have studied and parame- 
terized the transfer of o rbital to in ternal kinetic energy dur- 
ing galaxy encounters dRichstone 1975: Dekel et al. 198^; 
lAguilar fc Whitelll985h . lBinnev fc Tremainel (| 19871 ) summa- 
rize these studies and present formulas for impulse approxi- 
matio ns valid in the tidal and radial cases. lAguilar fc Whitg 
( 1985) find that the tidal approximation breaks down at ap- 
proximately 5i?e, but that a smooth interpolation between 
the two cases gives an approximate fit to actual energy trans- 
fers: 



REFERENCES 

Aguilar L. A., White S. D. M., 1985, ApJ, 295, 374 
Barnes J. E., 2002, MNRAS, 333, 481 



Tidal case 



Radial case 



AE 



2G^MiMir^ 

pcri peri 



AE - 



(16) 
(17) 
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Ml and M2 are the masses of the perturbed and perturbing 



Table 4. Definitions of the inputs used in the model. 



distance and velocity. For the tidal case, is the mean- 
square radius of the perturbed galaxy. A Plummer model, 
$ = —GM I \J r'^ + a^, is assumed for the radial approxima- 
tion. 

We measure the change in internal kinetic energy of 
each progenitor during the first close encounter and de- 
fine this as our "impulse." In order to compare with the 
above approximations we assume that the galaxies are point 
masses and use values of -Rperi and Vperi calculated at the 
beginning of the simulations. The impact parameters sim- 
ulated in our suite all fall well within 5-Rdm, as might be 
expected for merger orbits. Consequently, the tidal approx- 
imation performs quite poorly at predicting the impulse. 
Furthermore, any impulse approximation assumes that the 
dynamical time of the perturbed galaxy is much less than the 
time of the encounter. However, if two galaxies are going to 
merge, then the encounter velocity is typically of the same 
order as the internal velocity of the larger galaxy. There- 
fore, it is not clear that either impulse approximation above 
would apply. We compared the radial approximation, and 
several similar functions that included a dependence on im- 
pact parameter, to the measured impulses and found that 
they did not provide satisfactory fits. 

However, we do find that the measured impulse, during 
an encounter with orbital parameters within our range of 
simulated values, is well-fit by the following formula: 



'K?cri(-Rpcri + B -Rl.tot -Rpcii + C , 



(18) 



where Mi, tot and M2,tot are the total mass of the perturbed 
and perturbing galaxies respectively, baryons plus dark mat- 
ter, and -Ri^tot is the total half mass radius of the perturbed 
galaxy. The parameters are set by best fit to the simulation 
resuhs with A = 1.6, B = 1.0, and C = 0.006. 

For a radial orbit, this formula vaguely resembles the 
radial approximation. However, the mass dependence of the 
formula is different. In the radial impulse approximation 
A.E oc MiM|, whereas in our formula AE oc M1M2. This 
suggests that the impulse approximation is breaking down, 
that we cannot simply assume that the potential of the per- 
turbed galaxy is constant during the encounter. At inter- 
mediate impact parameters, where -Rpcri < -Ri.tot, the fit- 
ting formula falls ofi^ as l/i?pori. At larger impact param- 
eters, Rperi > Rdm, it falls off as l/Rp^ri- The parameter 
C sets the cutoff radius at which the impulse ceases to in- 
crease with decreasing impact parameter and approaches 
the constant radial case. The fit of C is somewhat tenu- 
ous since our simulations do not actually probe the range 
of Rperi < CRdm- However, these cases are so radial that 
the probability of such encounters is low, and the eff'ects of 
small changes in C are likely to be insignificant. Our suite 



probes a range of roughly 0.1i?i,tot < Rpai < 2i?i, 



and 



2crp 



< Vpo 



is the initial velocity 



dispersion of the larger progenitor. The formula is plotted 
against impulses measured from our simulations in Figure 
[3l This plot includes both equal and unequal mass mergers. 
The fit is valid over a range of progenitor mass distributions. 
Specifically, the structure of the fitted galaxies varies widely 
with a variety of bulge fractions and dark matter concentra- 
tions. 



Name 


Definition 


Inputs 


Mtot.l, Aftot,2 


The total masses (baryonic plus 




dark) of galaxies 1 and 2 respec- 




tively. 


Rtot 1 , Rtot 2 


The three-dimensional half-mass 




radii of the total masses (baryonic 




plus dark) of galaxies 1 and 2 re- 




spectively. 


^pcri ) ^pcri 


The theoretical pericentric velocity 




and distance of the first encounter, 




defined when the galaxies are sep- 




arated by ~ Rvir- 


K^tot.li ^tot,2 


The total kinetic energies 




of the galaxy/halo systems 




of galaxy 1 and 2 respec- 




tively. In our simulations 




/Ctot.i ^ 0.35 GMl^JRtot,! 




with very little scatter. 


/g.l' /g.2 


The gas fractions of galaxies 1 




and 2, defined as (gas mass) / (total 




mass). 


Ri,R2 


The stellar three-dimensional half- 




mass radii of galaxies 1 and 2. 


Mi,M2 


The stellar masses of galaxies 1 and 




2. 



A^dm.li ^^dm,2 The dark matter mass inside 1/2 
the stellar three-dimensional half- 
mass radius of galaxies 1 and 2 re- 
spectively. 

Afstars.l , Mstars,2 The stellar mass inside 1/2 the 
stellar three-dimensional half-mass 
radius of galaxies 1 and 2 respec- 
tively. 



8 APPENDIX B: SUMMARY OF THE MODEL 

The model presented requires the integration of formulas 
and ideas found throughout the paper. Hence, for prag- 
matic purposes, we present a brief summary of the model 
for the reader who wishes to implement it as a recipe within 
a semi-analytic model of galaxy formation. In Figure [5] we 
illustrate the outline of the model. The relevant equations 
are included, except for the definitions of the energy terms in 
the conservation equation. For these, we refer the reader to 
the list in i]3.3l Inputs, parameters, and outputs are defined 
in Tables H [S] and H 

The general flow of the model is as follows: 

1. Use orbits and masses to calculate impulse. 

2. Use impulse and gas fraction to calculate new stars 
formed in the burst. 

3. Use mass of new stars and properties of progenitors to 
calculate radius. 

4. Use initial dark matter fraction and final mass of stars 
to calculate remnant central dark matter fraction. 



Predicting the Properties of Merger Remnants 



Table 5. Definitions of the parameters nsed in tlie moclcl. 



Name/Value 


Description 


Parameters 


A = 1.6 
B = 1.0 
C = 0.006 


Parameters fit to match the 
impulse model to the simulations. 
See Appendix A. 


C'new — 0.3 


Proportionality constant in star 
formation ccjuation. Determines 
mass (jf iLvw slais. 


Cint = 0.5 


Structural constant which sets rel- 
ative weighting of internal energy. 


Crad = 1.0 


Constant which sets relative 
weighting of radiated energy. 


Csig = 0.30 


Structural constant which sots pro- 
portionality of (T^ to M/R in rem- 
nants. 


Cstars = 0.35 


Sets fraction of stars within 1/2 of 
the half-mass radius of the rem- 
nants. 



Table 6. Definitions of the outputs used in the model. 



Name 


Definition 


Outputs 




Mass of stars in the remnant (old 
+ burst). 


Ri 


Stellar three-dimensional half- 
mass radius of remnant. 


/dm.f 


The (lark matter fraction inside 
i/2 nf h'l 


(Tf 


The stellar velocity dispersion of 
the remnant inside the projected 
half-mass radius. 



5. Use central dark matter fraction and radius to calculate 
velocity dispersion. 
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Figure 9. Summary of the inputs (ovals), outputs (white rectangle), and equations (shaded rectangles) needed to implement the merger 
model. 



